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Abstract. The post-Minkowskian limit and gravitational wave solutions for general fourth- 
order gravity theories are discussed. Specifically, we consider a Lagrangian with a generic 
function of curvature invariants f(R, R a pR a ^, Rafiks R a ^ lS )- It is well known that when 
dealing with General Relativity such an approach provides massless spin-two waves as prop¬ 
agating degree of freedom of the gravitational field while this theory implies other additional 
propagating modes in the gravity spectra. We show that, in general, fourth order gravity, 
besides the standard massless graviton is characterized by two further massive modes with a 
finite-distance interaction. We find out the most general gravitational wave solutions in terms 
of Green functions in vacuum and in presence of matter sources. If an electromagnetic source 
is chosen, only the modes induced by R Q pR ai9 are present, otherwise, for any f(R) gravity 
model, we have the complete analogy with tensor modes of General Relativity. Polarizations 
and helicity states are classified in the hypothesis of plane wave. 
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1 Introduction 

Identifying the correct theory of gravity is a crucial issue of modern physics due to the fact 
that General Relativity, in its standard formulation, presents shortcomings at ultraviolet and 
infrared limits. For the first issue, we need a theory that should deal with gravity under 
the same standard of the other fundamental interactions (Quantum Gravity) [1]. In the 
other case, modifications of gravity are required to deal with the vast phenomenology coming 
from astrophysics and cosmology, generally addressed as dark matter and dark energy issues 
[2, 3], This matter is rather controversial due to the fact that the ambiguity comes out from 
the fact that phenomenology could be explained successfully both considering new material 
ingredients (dark matter particles addressing the problem of structure formation and light 
scalar fields giving rise to the acceleration of the Hubble fluid) or modifying gravity that, at 
scales larger than Solar System, could behave in different way with respect to the weak field 
limit related to General Relativity [4-8]. Furthermore, there are very few investigations and 
experimental constraints probing the gravitational field in very strong regimes. Several times, 
extrapolations of General Relativity are simply assumed without considering corrections and 
alternatives that could strongly affect theoretical and experimental results. 

With this situation in mind, it is urgent to find out some experimentum crucis or some 
test bed capable of discriminating among concurring gravitational theories [9-12], that, in 
any case, should reproduce the well-founded theoretical and experimental results of General 
Relativity. At astrophysical level, discriminations could come from anomalous stellar systems 
whose structures and parameters do not find room in the constraints and limits imposed by 
General Relativity. For example, extremely massive neutron stars, magnetars or compact 
objects like quark stars could be independent signatures for modified theories of gravity 
considered as extensions of General Relativity in the strong field regime [13, 14]. 

Besides, discrimination could happen in the realm of gravitational wave physics. This 
sector of physics, practically unexplored from the point of view of modified gravity, deserves 
a lot of attention since the large part of efforts has been devoted to the study of gravitational 
radiation in the realm of General Relativity discarding the fact that modified gravity presents 
a huge amount of new phenomenology and features. For example, only General Relativity 
strictly forecasts massless gravitons with two polarizations. In general, modified gravity 
and, in particular Extended Gravity, allows also massive and ghost modes and then further 
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polarizations [15, 16]. Specifically, several authors refer this issue to the Fierz-Pauli linearized 
analysis of massive gravity that leads to the so called van Dam-Veltman-Zakharov (vDVZ) 
discontinuity. In such a case, we are in presence of the Boulware-Deser ghosts. Such anomalies 
can be cured through the Vainshtein mechanism. A detailed review of these arguments can 
be found in [17]. 

However, the possibility of these massive modes are studied from a theoretical point of 
view but practically ignored from the experimental point of view due to the enormous difficul¬ 
ties related to the detection of gravitational waves. However, the forthcoming experimental 
facilities like VIRGO (after the Virgo supercluster of galaxies) [18], LIGO (Laser Interferom¬ 
eter Gravitational-Wave Observatory [19]) collaboration, LISA (aser Interferometer Space 
Antenna [20]) etc. could be suitable, in principle, for detecting these further modes. 

In this paper, we propose a systematic study of gravitational wave solutions in theories 
where generic functions of curvature invariants are considered generalizing the first partial 
outcome in the only f(R) framework [21], These are a straightforward generalization of 
f(R) gravity where the degrees of freedom, related to the curvature invariants, is considered. 
However, we have to say that we are not considering DR and similar terms where derivative 
of curvature invariants appear. We are also not considering the parity-odd Chern-Simons 
invariant [22] that enters at the same order in curvatures and derivatives. Here, we are 
taking into account only fourth order theories of gravity where derivatives of metric tensor 
appear up to the fourth ones. It is interesting to see that relaxing the hypothesis that 
gravitational interaction is derived only from the Hilbert-Einstein action, linear in the Ricci 
curvature scalar R, further gravitational modes, polarizations and helicity states come out. 
This new features are directly derived from the post-Minkowskian limit of the theory and 
points out a new rich phenomenology that deserves investigation in view of possible future 
detection of gravitational waves. The classification of these modes, implying new massive 
and polarization states, is the same of Wigner’s little group E( 2) [23, 24] earlier explored 
in literature by Eardley et al. [25]. Besides, f(R) theories of gravity have been very well 
explored in this sense. For example, it is well known that they are equivalent to scalar-tensor 
theories of gravity following the Bergmann-Wagoner formulation [26, 27] which generalizes 
the Brans-Dicke theory. Furthermore, kinematics of gravitational waves in f(R) gravity 
have been extensively studied by Berry and Gair which confronted gravitational radiation 
with Solar System tests [28]. It is important to stress that being f(R) gravity equivalent 
to a particular class of scalar-tensor gravity theories, the gravitational wave content has 
already been characterized in both representations [29]. Also a comparison between the two 
frameworks (the Einstein and Jordan frame) has been analyzed in the weak field limit [30, 31]. 
Finally, gravitational wave kinematics and dynamics in strong field regime and around the 
Minkowski background have been investigated in detail [32], 

This paper is organized as follows. In Section 2 we report briefly the field equations of 
fourth order gravity. In Section 3, we discuss the post-Minkowskian limit and the linearized 
field equations while, in Section 4, the gravitational wave solutions are reported. Section 5 is 
devoted to the discussion of all possible polarization and helicity states of the wave solutions. 
Conclusions are reported in Section 6. 

2 The field equations of Fourth Order Gravity 

The most general class of gravitational theories involving curvature invariants in four dimen¬ 
sions is given by the action 


- 2 - 


( 2 . 1 ) 


A 



f{X,Y, Z) + XL 


where / is an unspecified function of curvature invariants X = R, Y = R a pR a P , and 
Z = Ra^sR 01 ^ ■ The term L m is the minimally coupled ordinary matter contribution. In 
the metric approach, the field equations are obtained by varying ( 2 . 1 ) with respect to g^. 
We get the fourth-order differential equations 


HfiU = fxR^u - x9\iv - fx-!xv + g^Ofx + 2 fyR^Rav ~ 2 [fYR a {^\-v)a + □ [fyR^u] 

( 2 . 2 ) 

+ [fyR a ^ afi g^ + 2 fzR^Ru^ 1 - 4[/zR M a/3 J ;a/3 = XT fW 

where T^ v = — j s the the energy-momentum tensor of matter, fx = §x' 

fy = §f- fz = jyz-, □ = -a' a ■, and X = 87 rG 1 . The conventions for Ricci’s tensor is 

R^u = R a fiav and for the Riemann tensor is = T^ +_ The affinities are the usual 

Christoffel’s symbols of the metric: T^ = ^g 11 ' 7 (g a cr,p +g/3a,a — gap,a)- The adopted signature 

is (-|-) (we follow the conventions by Landau and Lifshitz [33]). The trace of field Eqs. 

( 2 . 2 ) is the following 


H = fxX + 2 fyY + 2 f z Z - 2/ + D[3f x + fyX\ + 2 [{fy + 2 f z )R a % a p = XT (2.3) 

where T = T a a is the trace of energy-momentum tensor and H = H a a . 

Some authors considered a linear Lagrangian containing not only X , Y and Z but also 
the first power of curvature invariants OR and R al3 - a /3- Such a choice is justified because all 
curvature invariants have the same dimension (L~ 2 ) [43]. Furthermore, this dependence on 
the two last invariants is only formal, since from the contracted Bianchi identity (2 R a ^- a g — 
□R = 0) we have only one independent invariant. In any linear theory of gravity (the 
function / is linear) the terms DR and R a/3 ]a /3 give us no contribution to the field equations, 
because they are four-divergences. However if we consider a function of OR or R a ^- a /3 by 
varying the action, we still have four-divergences but we would have the contributions of sixth 
order differential terms. As said in the Introduction, in this paper, we consider only fourth 
order differential field equations. This means that Action (2.1) is a fourth-order theory and 
no higher-order terms in derivatives are taken into account. However, as we will see below, 
one needs only two of the three curvature invariants, due to the Gauss-Bonnet topological 
invariant which fixes a constraint among the curvature terms. 

3 The Post Minkowskian Limit 

Any theory of gravity has to be discussed in the weak field limit approximation. This 
prescription is needed to test if the given theory is consistent with the well-established 
Newtonian theory and with the Special Relativity as soon as the the gravitational field 
is weak or is almost null. Both requirements are fulfilled by General Relativity and then 

1 Here we use the convention c = 1. 


- 3 - 






they can be considered two possible paradigms to confront a given theory, at least in the 
weak field limit, with the General Relativity itself. The Newtonian limit of f(R )-gravity 
and /(R, R a pR a ^, Ra/^R^^-gravity can be investigated always remaining in the Jor¬ 
dan frame [44-49] while a preliminary study of the post-Minkoskian limit for the /(R) 
gravity is provided in Ref. [50]. Here we want to derive the post-Minkowskian limit of 
/(R, R a ^R a ^, R a p 1 sR a ^ & ) gravity, that is for generic fourth-order theory of gravity, with 
the aim to investigate the gravitational radiation. 

The post-Minkowskian limit of any theory of gravity arises when the regime of small 
field is considered without any prescription on the propagation of the field. This case has to 
be clearly distinguished with respect to the Newtonian limit which, differently, requires both 
the small velocity and the weak field approximations. Often, in literature, such a distinction 
is not clearly remarked and several cases of pathological analysis can be accounted. The 
post-Minkowskian limit of General Relativity gives rise to massless gravitational waves and 
reproduces the Special Relativity. An analogous study can be pursued considering, instead 
of the Hilbert-Einstein Lagrangian, linear in the Ricci scalar R, a most general function / of 
curvature invariants. 

However, working with post-Minkowskian limit, one has to be extremely careful with 
the issue whether or not there are ghost modes that could render the theory meaningless. A 
standard way to deal with ghost modes is to consider any alternative theory as an effective 
field theory close to General Relativity, since it is well known that General Relativity has 
no ghosts. Ri fact, effective field theory has a restricted regime of validity, and the ghost 
modes only appear once one leaves such a regime. See, for examples [34, 35] where such a 
problem is discussed for the Chern-Simons and /(R) theories respectively. Essentially, the 
approach consists in decoupling the modes in the weak limit. In [32], it was shown that in 
the decoupling limit, and approaching the asymptotically flat region far from any matter 
sources, a very generic class of theories with higher curvature invariants leads exactly to the 
same gravitational wave modes as in General Relativity (see for example [36] for the case of 
f(T) teleparallel gravity). 

Another question is related the Cauchy problem of the theory, and in particular with 
its well-formulation and well-posed initial value formulation. For /(R) gravity, this problem 
has been highly debated and discussed [37-40]. In any case, considering theories with higher 
derivatives, the Cauchy problem has to be treated with care because the related Hamiltoni¬ 
ans could suffer, in general, with the Ostrogradski instability (see [41] for a detailed study). 
In particular, the existence of a well-posed initial value problem for higher-order theories 
of gravity involves the introduction of auxiliary degrees of freedom to reduce the derivative 
order. In some sense, the mechanism is the same acting in conformal transformations where 
the further degrees of freedom in the Jordan frame are disentangled in scalar fields minimally 
coupled to gravity in the Einstein frame [16]. This fact allows to evaluate how many prop¬ 
agating dynamical degrees of freedom are present in the theory. Alternatively, treating the 
theory as an effective field theory can contribute to address the question of the existence of 
a well-posed initial value formulation [42]. 

In order to perform the post-Minkowskian limit of field equations, one has to perturb 
Eqs. (2.2) on the Minkowski background rIn such a case, we obtain 


9fj.u — Vnv T h[n/ 

with hfj, v small (OQi^ 2 ) <C 1). Then the curvature invariants X, Y, Z become 
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(3.1) 


(3.2) 


'X ~ x« + e>(v 2 ) 

y ~ y^ 2 ) + o(v 3 ) 

~ z( 2 ) + o(v 3 ) 


and the function / can be developed as 


f(x, Y, Z) ~ /(0) + /*(0)X« + i/^ojX* 1 ) 2 + f Y (0)y ( 2) + /z(0)z(2) + o ( ^3 )( 3_ 3) 


Analogous relations for partial derivatives of / are obtained. From lowest order of field Eqs. 
(2.2) and (2.3) we have the condition /(0) = 0, while at 0(1) - order, we have 2 


fx(0)R$ + [/y(0) + 4/z(0)]D,ji2^ - ixMxWri,*,, + [fxx( 0) + ^]7 lflu n v X^ 

-/xx(0)I (1 V - [/r(0) + 4/ z (0)]i? (1) V a - /y(0)i? (1) V« = (3-4) 

-fx( 0)lW + [3 /x_y(0) + 2 /y (0) + 2/ z (0)]D r? X( 1 ) = 


where is the d’Alembert operator in the flat space. T^ u is fixed at zero-order in (3.4) 
since, in this perturbation scheme, the first order on Minkowski space has to be connected 
with the zero order, with respect to h^ u , of the standard matter energy momentum tensor. 
This means that T IIU is indipendent of h^ u , and satisfies the standard conservation conditions 
= 0. By introducing the quantities 


fx( Q) 


3fxx (0)+2/y (0)+'2 fz (0) 


(3.5) 


Sx{ Q) 


/r(0)+4/ z (0) 


we get two differential equations for curvature invariant A^ 1 ) and Ricci tensor Rf^ 3 



X« = m 2 2 XT$ 


(3.6) 


(□,, + mi 2 )X« = -mi 2 AT(°) 


We note that in the case of /(X)-theory we obtain only a massive mode (with mass mi) 
of Ricci scalar (X). In fact if f Y = fz = 0 from the mass definition (3.5) mi 2 —> 
(—3/xx(0)) _1 and m 2 2 —> oo we recover the equations of /(X)-gravity [50] 


2 We are using the properties: 2R a /3’ a P — OR = 0 and 

3 We set fx = 1 i-e. G -¥ fx(0)G. 
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(3.7) 


R 


(i) xw„ , dlu-n^n wn _ 

pu 2 'W “T 3mi 2 ^ — 


XT, 


(o) 


/it/ 


(□,, + mi 2 )xW = -rrn 2 XT^ 

while f(X, Y, Z)-theory we have an additional massive propagation (with mass m 2 ) of Ricci 
tensor. Finally in the case of / —> X also m\ 2 —> 00 and we recover the General Relativity. 

A first consideration regarding the masses (3.5) induced by f(X, Y, Z) -gravity is nec¬ 
essary at this point. The second mass m 2 is originated by the presence, in the Lagrangian, 
of Ricci and Riemann tensor square, but also a theory containing only Ricci tensor square 
gives rise to the same outcome. Obviously the same is valid also with the Riemann tensor 
square alone. Then such a modification of theory enables a massive propagation of Ricci 
Tensor and, as it is well known in the literature, a substitution of the Ricci scalar with any 
function of the Ricci scalar enables a massive propagation of Ricci scalar. We can conclude 
that a Lagrangian containing any function of only Ricci scalar and Ricci tensor square is not 
restrictive. This result is coming from the Gauss - Bonnet topological invariant Ggb defined 
as Ggb = X 2 — 4 Y + Z [51]. In fact by applying the variation with respect to the metric 
tensor g /iu in 4 dimensions to the quantity f d A Xy/—g Ggb = 0, we have 


<5 j d 4 xV^G GB = J dSy^giH^ - + H^r = j d*xV=9 H° B 5g>"( 3.8) 


= ^fr ,In four dimen- 


whprp H x2 — 1 s (V~9 x2 ) ttY _ 1 o T , f l f[Z _ U 

wneie Sg ^„ , n gl/ - Sg ^„ anu gg ^ v 

sions we have that is identically zero. For this reason the variation of the Gauss-Bonnet 
invariant generates a dimension-dependent identity, that is 


- 477^ + H% = 0 (3.9) 

By substituting condition (3.9) in Eqs. (2.2), at post-Minkowskian level, we find the same 
Eqs. (3.6) with a redefinition of masses (3.5). In the weak field limit approximation, we can 
consider as Lagrangian in the action (2.1), the quantity [48] 


f(X, Y, Z) = aX + bX 2 + cY 

Then the masses (3.5) become 


(3.10) 


m l = -2(3f+7) 

(3.11) 

m 2 2 = “ 

These masses have real values if the conditions a > 0, b < 0 and 0 < c < — 36 hold. This 
fact is extremely important in order to establishing a no-ghost constraint on such a theory. 
For gravity theories explicitly containing the Gauss-Bonnet term in the post-Newtonian limit 
see [52], 
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4 Gravitational wave solutions 


Once developed the post-Minkowskian limit, one can search for gravitational wave solutions. 
The general solution of field equations (3.6) is given by 


h[xu — T h/iv (4.1) 

where f)^ is the (homogeneous) solution in the vacuum and h Mi , is the (particular) one in the 
matter. First we try to find the solution From the second line of (3.6), by introducing 
the Green function Gkg,i(x,x') of Klein-Gordon field defined as follows 

(□jj + mi 2 )QKG,i{ x i x ') = 5 <yi \x-x'), with i = 1,2 (4.2) 

where — x') is the Dirac delta function in four dimensions, we find 


X (x) = ~ mi 2 X / d 4 x' g K G,i(x,x')T(°\x') 


(4.3) 


where x = x p = (t, x) = (t, x 1 , x 2 , x s ). The first line of (3.6) can be recast as follows 


(□ ri + m 2 2 )R$ = X + W ^ 


L |w 


6 


(m ‘ 2 ~ [ d * x - 


d 2 - —r, 

u fiu 9 '//•«' 


(4.4) 


Qkg,i{x,x')T ( 0 ) {x') 


so the solution for the Ricci tensor is obtained 


x j d 4 x' Gkg,2(x, x') m 2 2 T < p{x') - mi + 2m2 r]^ u T^\x') 


(4.5) 

—- I dV d A x" Gkg,2{x,x') d 2 , v , - Qkg,i(x' ,x")T^\x") 


The Ricci tensor, in terms of the metric (3.1), is given by 


R\iv ~ h \n,v)c 2 ^’^ ( 4 - 6 ) 

where h = h a a . Since we can use the harmonic gauge condition g pr7 T a pa = 0, we set 

hfj,a a — 1/2 h tfl = 0, then the Ricci tensor becomes i?)// = — ^\3 V h^. The solution of 
Eq.(4.5) is 
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h/ji/ — 


2(mi 2 — m 2 2 )X 


J (fix' (fix" (fix'" 9 gr(x, x') 9kg,2{x', x") 


r ) 2 

°H"v" 


m i 


-Vi* 


9kg,i(x",x'")T^\x"') 


(4.7) 


-2X J (fix' (fix" 9gr(x,x')9kg,2(x',x") 
where 9gr(x,x') is the Green function defined as 




O v 9gr(x,x') = 5 (4 \x-x') (4.8) 

Considering the expressions of the Green functions 9gr{x,x ') and 9kg, i,2(x,x') in terms of 
plane waves it is possible to rewrite the solution (4.7) as follows 


hnv(x) = X f d 


4 x' 


Z 2 {x,x') T${x!) + ( Z(x, x') r]^ + Z t 


■'flu 


(x,x')j T(°\x' 


h ( rc ) = X f (fix' Z\{x, x') T ( -°\x') 

where 


(4.9) 


Z s (x,x') 


2 m s 2 (-l) 1+s 


r d 4 k e >k(x-x’) 

J ( 2n ) 4 k 2 (k 2 —m s 2 ) 


Z(x, x') 


f d 4 k (m. 1 2 +2m 2 2 )k 2 -3m. 1 2 m.2 2 j k(x-x') 

J (27 r) 4 3 k 2 (k 2 —mi 2 )(k 2 —m 2 2 ) 


(4.10) 


Z^ u (x,x') 


C d 4 k 2(iiir-inZ)k ll k,, j k(x-x') 
J (2tt) 4 3 k 2 (k 2 —mi —m 2 2 ) 


k = = ( w , k ) = ( a ;, k 1 , k 2 , k 3 ), kx = k a x a = ut — k • x , k 2 = k a k a = a ; 2 — | k | 2 

and s = 1, 2. We note that in the case / —> X (from the mass definitions (3.5) we have 
mi, m 2 oo) we find Z s (k) —> 2(—1 ) s k~ 2 , Z[k) —k 2 , Z fiu (k) —> 0 and the solutions 
(4.9) become those of General Relativity: 


h^ix) = -2X J cfix'9 G R(x,x') S$(x'), (4.11) 

where S^J = T$ - rj^T^/2. 

Finally we can recast the propagators (4.10) in terms of the Green functions 9gr{x,x') 
and 9kg,s(x, x') and we can see immediately the propagation of the massless interaction and 
of two massive interactions. We get 


Z 8 (x,x') : 
Z(x,x') = 
Z tiu (x,x') 


2 (— 1 ) 


1+S 


9gr{x , x') - 9 kg, s (x, x') 


9gr{x,x') - \ 9kg,i{x,x') - \9kg,2{x,x') 


2 d 2 
3V 


9gr(x , x') + g ^Ax,x') _ g KG , 2 (x,x') 

m± Z m2 Z ^ ' V ’ ' 777,1 z 777,0^ 


777-1 z 


m 2 z 


(4.12) 
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It is important to stress that the limit m —> 0 could become problematic for the above 
prapagators since the theory appears strongly coupled in this limit. However the source term 
can be assumed decoupled as soon as m —> 0. In such a case the problem is solved. See [17] 
for a discussion. 

A very interesting feature is achieved when we consider a traceless source with T 1°) = 0. 
In this case, solutions (4.9) are unaffected by the functions of the Ricci scalar but only by 
the presence of the Ricci tensor square terms into the interaction Lagrangian. Moreover, 
from the second line of (4.9), we have also that the gravitational waves must be traceless. A 
prototype of traceless source is the electromagnetic field where we have 


h^(x) = X J d 4 x'Z 2 (x,x')T$’ em (x') with h em (x) = 0. (4.13) 

Finally, the solution of the first field Eqs. (3.6) in the vacuum (T ^ = 0) can be derived. By 
using again the hypothesis of harmonic gauge and the principle of plane wave superposition, 
we get 




2 2 
777, i — 777-2 

3 m 2 2 


d 4 x' Z- 2 (x, x') 


r/5 2 

mi 2 



-2 


d 4 x'g GR (x, x')R { ^ hs) (x') + h^^ hs) 


(4.14) 


where X^, h^(hs) are respectively the homogeneous solutions of the Klein-Gordon 

equation for the Ricci scalar and tensor and the solution of wave equation for the metric. The 
homogeneous solutions are chosen in such a way that they satisfy the boundary conditions 
and the gauge harmonic condition — 1/21)’^ = 0. Also here, the limit m —> 0 becomes 
problematic since it is strongly coupled. As above, the remedy is that the sources decouples 
in this limit [17]. 

An important remark is necessary at this point. In the above calculations, we are using 
the harmonic gauge condition that allows to transform admissible field configurations if one 
considers sufficiently small regions of spacetime. However, such a condition has to work 
over more than one gravitational wavelengths and the effects of background curvature on 
wave propagation have to be considered. A detailed discussion on harmonic gauge condition, 
shortwave approximation and propagation on curved background can be found in [53]. Here 
we adopted the same approach (see § 35.13, § 35.14 therein). 

Below we will consider polarizations and helicity states in vacuum starting from this 
result. 


5 Polarizations and helicity states in vacuum 

We can analyze the propagation in vacuum by performing the Fourier analysis of field Eqs. 
(3.6). In the Fourier space, we have 


k 2 < (m 2 - Ar) 


9 9 

777,2 —mi 

3?77i 2 


kfiku + 77 



mi 2 +2777,2 2 

6777,1 2 



= 0 
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A 2 (mi 2 — k 2 ) f) = 0 


(5.1) 










where f) /i; ,, f) are the Fourier representation of f)^, f). The gauge condition now becomes 
fyfur k a — 1/2 fl k^ = 0. The solutions are shown in Table 1. 


Case 

Choices 

Gauge condition 

A 

k 2 = 0 

any with f) = 0 

o 

II 

b 

b 

B 

k 2 = mi 2 

. Tj i_u/ ^ 

9/"' — YvX 2 3 

with f) / 0 

verified 

C 

k 2 = m-2 2 

any with f) = 0 

o 

II 

b 

b 

D 

k 2 = 0 

any with f) / 0 

i)»a k a - 1/2 jj k„ = 0 


Table 1. Classification of solutions in the vacuum of field Eqs. (5.1). Specifically, the case A is one 
of General Relativity 

; the case B is massive and pure trace; the case C is massive, trace-free and transverse and 
the case D is massless where the trace-reverse of 1) is transverse. 

The case A corresponds to the standard massless gravitational waves of General Rela¬ 
tivity. The solution B is the same of the pure /(X)-gravity, i.e. the choice k 2 = my reduces 
automatically the field equations of f(X, Y, Z)-gravity to those of f(X)- gravity. Cases B and 
C are very different. In fact the solution C is traceless while the solution B satisfies directly 
the gauge condition. 

The solution of case B for 1)^ is given by the following expression 


h B _ 


d A k 

J^T 4 


kfjkjj ti^ u 

mi 2 2 


1 f{k)e' kx = - 


2 


d 2 

mi 2 


fhw 

2 


d 2 

m\ 2 


b B (x) , 


d A k 

(2 W 


tj B (k) e 


j kx 


(5.2) 


where the trace of metric is a generic Klein-Gordon function with k = (aq, k), where uj\ = 
\/|k| 2 + mi 2 . We can write the general solution in terms of its Fourier modes which are 
plane waves, that is 


I= / |Xc(k) (5.3) 

where C(k) is the Fourier representation of the trace t) B . If we consider a propagating trace 
in the z-direction 4 z) = f)o e’ z ), where k = (uq, 0, 0, k z ) with uq 2 — k z 2 = mi 2 , 

the solution (5.2) is given by 

4 We set x 3 = z. 
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(5.4) 



f B J (ui t-k z z) 


to 
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/I+ iV 

2 ^ mi' 

0 

0 

_ un k z 

mi 2 


0 

1 

2 

0 

0 


0 

0 


Uii k z \ 
-T" 


mi' 

0 


J+ fc* 

2 "l" mi 2 


e ) {uit-k z z) 


where e^, is the polarization tensor. By a change of coordinates x p —> x l,J ‘ = x p + (, p (x) with 
0(C 2 )< 1 , we can transform the metric into a new metric f)' = — C^/x • Let 

us suppose that we choose (^(x) = )9 fJ 'e ]kx , the metric t) /w (x) becomes = e'^ w e )kx 

where e' jJV = e^ v + k lt 0 u + k v 6^. By performing a change of coordinates and choosing 6^ = 
= (“i^T “ 2 mi^’ °’ °> 2 ^ ~ 4 ^ 72 ) th e polarization tensor e% v in Eq. (5.4) becomes 



/o 0 0 0 \ 


/0 0 0 0\ 


/0 0 0 0\ 

f)0 

0 0 0 

f)o 

0 10 0 

k z 2 f)o 

0 0 0 0 

3 

0 0-^ 0 

6 

0 0 10 

6 oil 2 

0 0 0 0 


0 0 --2+i^J 


\o 0 0 1 ) 


\0 0 0 1/ 


In the case C, we have 



CV(k) e L^*- k ' x ) 


(5.6) 


where C fiU ( k) is the Fourier representation of the gravitational wave and U 2 = y |k| 2 + m 2 2 . 
Also in this case, by considering a propagating wave in the ^-direction z) = ^ i u ^t-k z z) 

where the polarization tensor e^ u satisfies the traceless condition r] pa e^ a = 0 and after also 
the harmonic gauge ejy k a = 0 (in the Fourier space), we find the solution 


/ 

£00 

£01 

£02 

-f e oo\ 


£01 

£11 

£12 

-f e oi 


£02 

£12 

m 2 2 e 

~~jTT e 00 - £ll 

“fc7 e 02 

v- 

" fcl e °° 

w 2 , 

~fcT e oi 

“fc7 e 02 

Wo 2 f 

k7 e 00 / 


(5.7) 


where eoo, £oij £ 02 j £ 11 , £12 are unspecified values. By performing also in this case a change of 
coordinates [15, 16] and choosing 9 M = 9^ = ( — the polarization 

tensor e ( jnJ in Eq.(5.7) becomes 
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/° 

0 

0 

0 

\ 

0 

£ll 

^12 

m 2 2 

k z 0 J 2 

£01 

0 

£12 

-fyreoo - £li 

m2 2 

k z oj 2 

£02 

0 

m 2 2 , 

m 2 2 c 

m 2 4 


\ u 

k z U2 

k z ui 2 C ^ z 

k z ' 2 U2 2 

£00/ 



/0 0 0 0\ 


/o 0 0 0\ 


0 10 0 


0 0 10 

£11 

0 0-10 

+ £12 

0 10 0 


\0 0 0 0^ 


\0 0 0 0/ 



/o 0 0 0\ 


/o 0 0 0\ 

m 2 2 £oo 

0 0 0 0 

m 2 4 e 0 o 

0 0 0 0 

k 2 

0 0 10 

k z 2 u 2 2 

0 0 0 0 


\0 0 0 0^ 


\0 0 0 1/ 


Finally for the case D we have 



/o 0 0 0\ 


/0 0 0 0\ 

m-2 2 £01 

0 0 0 1 

m 2 2 £ 02 

0 0 0 0 

k z U 2 

0 0 0 0 

k z U 2 

0 0 0 1 


\o 1 0 0^ 


\0 0 1 0/ 



A 

— £03 

£01 

£02 

£03 


£01 

£ll 

£12 

—eoi 


£02 

£12 

-£ll 

—£02 

V 

£03 

—£01 

— £02 

-e 


(5.9) 


where eoi, 602 , £o 3 ) en, £12 are unspecified values and e is the trace of polarization tensor e^ u . 
Here if we choose 0^ = 9^ = (, — yy, — yy, — Al^ +£ ) the polarization tensor in 
Eq.(5.9) becomes 


(0 

0 

0 

0 ^ 


/o 0 

0 

0\ 


/0 0 0 0\ 

0 

eii 

£12 

£01 + £13 


0 1 

0 

0 


0 0 10 

0 

£12 

— £ll 

£02 + £23 

= £11 

0 0 

-1 

0 

+ £12 

0 10 0 

\0 

eoi + £13 

£02 + £23 

0 ) 


^0 0 

0 

0/ 


\0 0 0 0/ 


(5.10) 


+(eoi + £ 13 ) 


/0 0 0 0\ 
0 0 0 1 
0 0 0 0 
\0 1 0 0/ 


+ (£02 + £ 23 ) 


/0 0 0 0\ 
0 0 0 0 
0 0 0 1 
\0 0 1 0/ 


However the last two polarizations in (5.10) are not physical and pratically we obtained the 
same outcome of General Relativity (case A). 

If we introduce an independent basis (see Table 2) for the polarizations (5.5) and (5.8), 
the general solution of field Eqs. (3.6) in vacuum is 


fl z ) = 


'Hi + U 2 


(x) 


j k z (t-z) 


+ HL 3 


=( 1 ) 


k z 2 C (S) 




J (u>it—k z z) 


+ 


^ 4 £ 


(+) 


+ + W 6 (V3 £ « - V2'M _ ^±^1 e P>) (5.11) 


T H.7 e^ 2 J + ^8 £ 


(3) 

/il/ 


J ( 1 u 2 t-k z z ) 
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where Hi, H 2 are arbitrary constants related to the propagation modes of gravitational 

Hi = eii, 


waves in General Relativity and the other constants are defined as Hi = — v 'fo° . 


H 5 = eia, H 6 = 


m 2 2 eoi ^ _ v / 2m 2 "(;o2 


k z ui2 


k,ui 2 



/o 0 0 0\ 


/o 0 0 0\ 


/o 0 0 0\ 

1 

0 10 0 

,(x) _ 1 

0 0 10 

f s ) _ 

0 0 0 0 

V2 

0 0-10 


0 10 0 


0 0 0 0 


\0 0 0 0 ) 


\0 0 0 oj 


\0 0 0 1/ 



/o 0 0 0\ 


/o 0 0 0\ 


/0 0 0 0\ 

1 

0 10 0 

,(2) _ 1 

0 0 0 1 

,(3) _ 1 

0 0 0 0 

V3 

0 0 10 

e ^ u - s/2 

0 0 0 0 

e » v - s/2 

0 0 0 1 


\o 0 0 1 ) 


\o 1 0 Oj 


\0 0 1 0/ 


Table 2. The basis of the polarizations. Each polarization satisfies the condition e pa e pa = 1. 

The different components of polarization tensor can be distinguished if we ask how e pu 
changes when the coordinate system undergoes a rotation of a given angle ip about the 2 -axis 
[54], This is a Lorentz transformation of the form 


/I 0 0 0\ 

= 0 cos»> sin, 0 (5.12) 

^ U — sm ip cos ip (1 

\0 0 0 1/ 

and since it leaves k^ invariant (7 Zf 7 k a = kf), the only effect is to transform einto 
e pi/ = TZ P P 7Z u a e pa . In the case B, the polarization tensor e'f is unchanged, then we can 
state that the helicity is null 5 * . In the case C, we have 


■lC _ e ±j 2 ip pc m 2 gQO 


c ± 


H-^ 5 — sm <p e 

k 


±j (95+71-/2) 


If 


5 ±]ip lf 


(5.13) 


where e± = ep =F j £12 and l± = ei 3 =F j £ 23 - The helicity is ± 1 for the states l±, while for e± 
we cannot identify the condition if = ^ ef. As such condition holds only if eoo 7 ^ 0 

(and He 7 ^ 0 in Eq.(5.11)). In other words, we have the possible helicity values 0, ±1, ±2. 
Clearly, all the polarization tensors discussed here can be referred to the Wigner little group 
E( 2) as discussed in [25]. 




5 Any plane wave ip transforming under a rotation of an angle 95 about the direction of propagation into 

= e'^^ip has helicity £. 
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6 Conclusions 


In this paper, we presented a complete study of post-Minkowskian limit and gravitational 
wave solutions of fourth-order gravity theories of gravity in four dimensional space-time. We 
considered a generic action constructed with curvature invariants derived from the Riemann 
tensor, that are the Ricci curvature scalar, the squared Ricci tensor and the squared Riemann 
tensor. Thanks to the Gauss-Bonnet topological invariant, it is possible to consider only two 
curvature invariants since the third is always related to the others by the Gauss-Bonnet 
constraint. In this sense, all the budget of degrees of freedom can be related to R and R a p- 

With respect to standard General Relativity, new features come out in the post-Minkowskian 
limit. First of all two further massive scalar modes emerge in relation to the non linearity 
in the Ricci scalar and tensor terms. This means that massive gravitons are a characteristic 
of this theories and their effective masses are strictly related to the form of the action /. 

As supposed by several authors, these particles could play a fundamental role for the dark 
matter issue that, in this case, could directly come from the gravitational part of cosmic 
dynamics [55-57]. 

Furthermore, the total number of polarizations is six and helicity can come into three 
distinct states. It is worth noticing that we have not to choose arbitrarily x and + polar¬ 
izations as in General Relativity but all possible polarizations naturally come out. This fact 
could be of great interest for the gravitational wave detection since running and forthcoming 
experiments could take advantage from this theoretical result and investigate new scenarios. 

In a forthcoming paper, a detailed study of sources compatible with these results will be 
pursued. 
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